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ON STATION HEIGHT DETERMINATION 
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ABSTRACT 

The effect of omitting the parallactic refraction correction for satellite 
optical observations in the determination of station coordinates is analyzed for 
a large satellite data distribution. A significant error effect is seen in station 
heights. A geodetic satellite data distribution of 23 close earth satellites, con- 
taining 30,000 optical observations obtained by 13 principal Baker-Nunn camera 
sites, is employed. This distribution was used in a preliminary Goddard Earth 
Model (GEM 1) for the determination of the gravity field of the Earth and geo- 
centric tracking station locations. 

The parallactic refraction correction is modeled as an error on the above 
satellite data and a least squares adjustment for station locations is obtained 
for each of the 13 Baker-Nunn sites. Results show an average station height 
shift of +8 meters with a dispersion of ±0.7 meters for individual sites. Station 
latitude and longitude shifts amounted to less than a meter. 

Similar results are obtained from a theoretical method employing a proba- 
bility distribution for the satellite optical observations. 
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EFFECT OF PARALLACTIC REFRACTION CORRECTION 
ON STATION HEIGHT DETERMINATION 


I. INTRODUCTION 

Parallactic refraction correction for optical observations is modeled as an 
error on a large geodetic satellite data distribution. The effect of this error on 
the determination of station coordinates is analyzed. A description of the paral- 
lactic refraction correction is presented in the introduction including a simpli- 
fied error effect on station height. 

Two methods of error analysis are applied. Method A employs an actual 
satellite data distribution on 23 close earth satellites for 13 stations. The data 
for this case is considered as a discrete sample. Method B employs a satellite 
probability distribution for observations at an idealized station. The distribution 
is assumed to represent the theoretical population for the discrete sample of 
Method A. 

Data distributions and analysis techniques are discussed in Sections II and 
IH. Results are presented in Section IV and a special case of interest for paral- 
lactic refraction correction is discussed in the appendix. 

1. Parallactic Refraction 


In the plate reduction of a photographic observation of an object, the observa- 
tion is adjusted for normal atmospheric refraction by treating it as a star image. 
For a close earth object such as the satellite, see Figure 1, this amounts to an 
over correction for refraction. In the figure the parallactic refraction, A R, ad- 
justs the normal refraction R so that the direction of the observation is pointed 
to the satellite. 

The satellite optical observations are given in terms of right ascension (a ff ) 
and declination ( S ) which define the direction of the satellite star image as 
shown in Figure 1. 

The satellite zenith, Z, is 


Z = ^ - A R. 


( 1 ) 


1 



Apparent Satellite 
Direction (Before 
Refraction Correction (R)) 


STATION 

ZENITH 



(a 0 ,5 0 ) — Right ascension (a 0 ) and declination (S 0 ), define the direction 
of satellite star image corrected for normal refraction R. 

Z 0 — Zenith angle corresponding to (a 0 , 5 0 ) . 

R — Normal optical atmospheric refraction. 

AR — Parallactic refraction. 


Figure 1. Parallactic Refraction 

The parallactic refraction correction used in the analysis as a modeling error 
is 


AR = 2.1 tan Z /p cos Z (radians) 


( 2 ) 


where p is the distance (meters) from the station to the satellite. A more 
complete representation accounting for variations in atmospheric conditions and 
short range distances is given in references 1 and 2. 

2. Station Height Adjustment 

If the station is constrainted from any horizontal adjustment, a simple re- 
lationship exists between parallactic refraction error and station height adjustment 


2 






as illustrated in Figure 2. The relationship is given in equation (3) and results 
in a formula, equation (4), of height adjustment (Ah) as a function of zenith angle 
only. The formula is shown to be applicable in Section II, 2.1 for the case of 
the theoretical satellite data distribution of Method B, because of symmetry in 
satellite geometry about the local vertical of the station. In the case of the dis- 
crete sample of Method A for an actual satellite data distribution, no assumptions 
are made about the distribution and all three components of the station coordinate 
are allowed to adjust. 


LOCAL 

VERTICAL 



Figure 2. Station Height Adjustment 


Using the law of sines for the triangle in Figure 2 and the fact that AR is 
small, 


A h = p A R/sin Z (meters) (3) 

where p is in meters andAR in radians. From A R in equation (2) above, equation 
(3) becomes 


Ah = 2.1/cos 2 Z (meters) 


( 4 ) 


Table 1 is presented to show values of Ah and AR as a function of zenith angle. 
For AR, p = 1000 kilometers is used for simplicity. It is seen that Ah increases 
with zenith angle to a value as large as 18 meters for zenith Z = 70°. 

The height relationship of equation (4) is used in Method B to estimate an 
expected height adjustment Ah over the theoretical satellite frequency 
distribution. 
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Table 1. Values of Ah and AR As a Function of Zenith 


z° 

Ah (meters) 

AR (seconds'of arc, p = 1000 km.) 

30 

2.8 

0.3 

40 

3.6 

0.5 

50 

5.1 

0.8 

60 

8.4 

1.5 

70 

18.0 

3.5 


Values of A R as seen in Table 1 are of the order of size of the accuracy (a-) 
of optical observations, namely 2" of arc. In view of this a separate random 
error from a normal population (a — 2") is applied to each of the observation 
components of right ascension and declination for the case of the discrete 
sample of Method A. 
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II. SATELLITE DATA DISTRIBUTIONS AND ANALYSIS - METHODS 
A AND B 

1. Method A - Actual Satellite Data Distribution (Discrete Sample) 


The parallactic refraction AR in equation (2) is modeled as an error on 
satellite optical observations of right ascension and declination ( a , S ). The 
effect of this error is analyzed on the determination of station coordinates for 
a large geodetic satellite .data distribution. The satellite data distribution con- 
sists of 30,000 optical observations (a, S ) taken by 13 Baker-Nunn camera sites 
on 23 close earth satellites, ranging in altitude from about 400 to 5000 kilometers. 
A least squares adjustment in station latitude, longitude, and height are deter- 
mined for each of the 13 stations. The mathematical equations for the model 
are developed in Section III. The satellite data distribution is presented in 
Table 2 and its distribution as a function of zenith angle is presented in Figure 3. 
Figure 3 is discussed in the next section where a comparison is made between 
the discrete and theoretical frequency distributions for zenith angle. 



Figure 3. Satellite Frequency Distribution in Zenith Angle 
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Table 2. Satellite Data Distribution for Method A 
(Baker-Nunn Data) 


Satellite 

No. of 
OBS 

Mid 

Altitude 

Perigee 

Height 

Apogee 

Height 

(Kilometers*) 

-AGENA-R 

640011 

477 



950 

ANNA-1 E 

620601 

1630 



1200 

BE-B 

640841 

126 


UBS* : ■ 

1100 

BE-C 

650321 

2286 

1100 


1300 

COURIER 

600131 

1447 

1050 


1200 

DI-C 

670111 

425 

950 

600 

1350 

DI-D 

670141 

2999 

1250 

600 

1900 

ECHO 

600092 

1014 

1600 

1500 

1700 

G EOS-1 

650891 

5812 

1700 

1100 

2300 

GEOS-2 

680021 

3660 

1350 

1100 

1600 

GRS 

630261 

174 

850 

400 

1300 

INJUN-1 

610162 

321 

950 

900 

1000 

MIDAS-4 

610281 

4451 


3500 

3750 

OGO-2 

650811 

363 

| ; • . | 

450 

1500 

OSCAR-7 

660051 

622 

1000 

850 

1200 

OVI-2 

650781 

315 

1950 

400 

3500 

SECOR-5 

650631 

124 

1800 

1150 

2450 

TELSTAR 

620291 

1620 

3300 

950 

5600 

TRANSIT 

610151 

728 

950 

900 

1000 

VAN-2R 

590012 

34 

1950 

550 

3300 

VAN-2S 

590011 

472 

1950 

550 

3300 

VAN-3S 

590071 

329 

2100 

500 

3700 

5 BN- 2 

630492 

178 

1100 

1050 

1150 


Heights rounded to nearest 50 kilometers 


The satellite data distribution in Table 2 was used in a preliminary Goddard 
Earth Model GEM 1, reference 3. A special case of interest in connection with 
parallactic refraction correction for the GEM 1 data distribution is presented in 
the appendix. 

2. Method B — Satellite Probability Distribution 


The satellite data distribution of Method A is regarded as a discrete sample 
from a theoretical population. 
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The probability distribution assumes the satellite population has equal area 
probability on any given sphere centered at the station, with a limiting condition 
on scheduling the satellite to within 75° in zenith angle. The probability 
distribution has complete symmetry about the local vertical of an idealized 
station. 

Consider the line of sight to the satellite at zenith Z and distance p as in 
Figure 2. Form an element of area dA at zenith Z in a ring of width dz on the 
sphere of radius p, then 


d A = 2 v p 2 sinZdZ 


Thus the probability (P) for the satellite to occur at zenith Z is 


d P = 


dA sin Z d Z 


I 


75 


sinZdZ 


The satellite frequency distribution is then 


( 5 ) 


(6) 


f (Z) - 


s in Z 

1 - cos 75° 


(?) 


The distribution f(Z) is plotted in Figure 3 for comparison with the histogram 
of zenith frequency for the discrete sample of Method A. The histogram is given 
in terms of percent frequency at 10° zenith intervals. In order to compare with 
the histogram, the theoretical frequency f(Z), which has units of percent per 
radian, is plotted with a scale factor of 10° in radian measure. The behavior of 
the sample case as given in the histogram is irregular beyond 70° and is plotted 
at 5° intervals for 70° to 80° to show a negligible frequency beyond 75°. The 
probability that the satellite will have a zenith angle as large as Z is equal to 
the area under the frequency plot from zero out to Z (after the above factor of 
10° in radian measure is inversely applied). Thus the integrated frequency 
distribution over the domain of Z for each plot will equal to unity. 

It is shown analytically in section 2.1 that use of equation (4), where Ah = 

2.1 /cos 2 Z, is applicable for Method B. This is so because of complete 
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sy mm etry in the satellite probability distribution about the local vertical at the 
station. The expected value of A h, over the frequency distribution in equation 
(7) is 


Ah 



A h f (Z) d Z 


(8) 


2.1 

cos 75° 


8. 1 meters. 


The following table is presented to show variation in A h for a 2° offset on 
each side of the limiting zenith angle of Z c = 75°. 


Table 3. Variation of Ah for 
Selected Zenith Limits 


Zc 

Ah meters 

73° 

7.2 

75° 

8.1 

77° 

9.3 


The table shows sensitivity to the zenith limit, but a value of 75° appears to 
serve the best viewing limit based upon the histogram in Figure 3. 

2.1 Simplified Height Relationship for Method B 

The simplified formula used in Method B for Ah = 2.1/cos Z resulted from 
the relationship shown in Figure 2, Section I, between the parallactic refraction 
A R and Ah. In this relationship it was assumed, as mentioned above, that no 
adjustment was necessary in the lateral direction of the station for the theoretical 
distribution of Method B. This condition is shown to be valid, although it may be 
intuitively clear. With the aid of Figure 2 consider a second point of satellite 
position placed symmetrically to the left of the station's local vertical. These 
two points of satellite position have equal probability of occurring in Method B. 
Through use of Figure 2, denote the vertical and lateral coordinates of the two 
satellite positions by (y, ±x) and denote the station coordinates by (v, %). From 
the relationship, 


tan Z = (x - l)/(y - v) 


(9) 
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the variation in Z gives respectively for the two points (y, ±x) of satellite 
position the two equations 


A Z = A R = 


s in Z 
~~P 


A v + 


cos Z M 

P 


(10) 


where AR = 2.1 tan Z / p cos Z from equation (2). 

Solving these two equations give A^ = 0 and the desired result for Ah , 
where A h = A v . Since the vertical plane in Figure 2 is arbitrary along with 
an arbitrary p and Z, the above result is valid for the entire distribution of 
Method B. 
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III. MATHEMATICAL MODEL FOR METHOD A 
1. Method A - GEM 1 Satellite Data Distribution 


The GEM 1 satellite data distribution, was described above in Section II. 
Table 2 above presented the distribution of data by satellites and Figure 3 gave 
the observation frequency distribution in zenith angle. 

For the error model, observation residuals for right ascension (a) and 
declination ( S) are formed to contain the error of parallactic refraction AR . 
The residuals include a random error from a normal distribution with standard 
deviation of 2 seconds of arc. These residuals are expressed in terms of a 
linear adjustment in station components. A least squares reduction is then 
performed on the satellite data for each of 13 stations. The satellite position 
relative to the stations was obtained from the satellite orbits reduced in the 
GEM 1 solution, reference 3. 

The observation residual equations in the form of orthogonal components 
are formally 


e 


l 


+ cos § A a - cos S (— — A h + A N + A (H) 

\3h B N BE / 

e xAS^liAh+^-iAN + ii-AE ( 12 > 

2 B h BN BE 


where 


A a. = a Q - a , transformed effect of A R 
AS = S Q - S , transformed effect of A R 
£ i> e 2 ~ the random errors described above 

Ah, AN, AE - orthogonal station component adjustments in metric length, respec- 
tively in the directions of the local vertical, North and East. 

The equations for transforming AR to the observation residuals for a, S and 
the equations for the partial derivatives may be developed with the aid of the 
geometry presented in Figure 4. The satellite/station orientation in the figure 
is presented on a unit sphere where arcs of great circles are employed. The 
angles defined in Figure 4 d enote the followin g quantities: 

Preceding page blank 
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Direction of Earth's 
Rotation Axis 



4 > - station latitude 

A z - satellite azimuth as given in the local horizon plane at the station 
and measured eastward from true north 

Z - satellite zenith 

Z 0 - Z + A R modeled observation of satellite with parallactic refraction 
error A R 
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a - right ascension of satellite measured eastward from Aries (7) 

S - declination of satellite measured northward from equatorial plane. 

t - angle between the vertical planes of the satellite and station and 
measured in the direction of the right ascension a of the satellite 

q - parallactic angle used to provide the orientation of AH and its pro- 
jections in the direction of increasing S and a . 

Thus from the definition of q as given above and as shown in Figure 4, the 
orthogonal observation error components due to A R are 

cos S (Oq - a) - A R sin q (13) 

A S = - A R cos q. (14) 

From the spherical triangle in the figure, the laws of sines and cosines give 


sin q = sin t cos 0/sin Z (15) 


cos q = (sin 0 - cos Z sin 8)/sin Z sin 8. (16) 

Denoting the station and satellite geocentric vectors by r and s respectively, 
then the satellite observing direction relative to the station is along the vector 


p = ' r - s. 


(17) 


The variation of (17) due to the effect of parallactic refraction error on station 
adjustment is 


Ajo= - As. 


(18) 


The orthogonal metric components of A p in the observation coordinate system 
of (o, a, 8 ) and As in the local coordinate system (V, N, E) defined above are 
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A p - (A p, p AS, p cos 8 A a) 

(19) 

A s = (A h, AN, A E) 

(20) 


Using Figure 4 the components of - As may be transformed to the com- 
ponents of A p through the following rotation matrices, 


A p 


Ah' 

p AS 

= - [*]■ [t] [- 4>] 

AN 

p cos 8 A a 


_A E_ 


(21) 


where the angles <p, t, and S are defined above. The matrix relation of interest 
in (21) for the observation equations (11) and (12) is given as follows: 


cos S A a 

1 

-St c <p c t 

s t s 0 


P 



A S 


cSsc£-sSctc<£ -stsS 

sScts0+c8c0 


where s - sin and c = cos for the trigonometric functio ns . 


Ah 

A E 

A N_ 
(22) 


Equation (22) provides the observation residual equations for the least 
squares adjustment of the station coordinates, for height Ah, for A N in the direct 
tion of latitude, and for A E in the direction of longitude. As indicated in equations 
(11) and (12) random errors el and e2 are added to the above residuals in the 
reduction (since the parallactic refraction error is the order of accuracy of the 
observations). 


Results are obtained for each of the 13 Baker-Nunn stations and presented 
in Section IV. 
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IV. RESULTS 

Error anal ysis results are presented for the effect of parallactic refraction 
error on station coordinates. Results for 13 Baker-Nunn stations are given in 
Table 4 using the technique of Method A, where the GEM 1 satellite data distri- 
bution was employed. Station height adjustments show an average or mean shift 
of 8.1 meters with an rms of 0.7 meters about the mean value. Horizontal com- 
ponents show very little adjustment. In Table 4 the mean value for the station 
height adjustment is the same as that for the idealized result of Method B, 
where a limiting value of zenith angle of 75° was employed. 

The small horizontal station adjustments of Method A reflect a good 
geometrical distribution of data for each of the 13 Baker-Nunn sites in the GEM 1 
data sample. In the idealized distribution of Method B it was shown in Section 
II, 2.1 that no adjustment results in the horizontal station components from the 
error effect of parallactic refraction correction. 
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Table 4. Station Adjustments for Method A and Method B 


Method A 


Station 

0° 

(Latitude) 

A h(m) 

AE(m) 

A N(m) 

Name 

(Height) 

(Longitude) 

(Latitude) 

1. ORGAN 

32 

8.0 

-0.1 

0.8 

2. OLFAN 

-26 

8.2 


0.0 

3. WOO MR 

-31 

10.2 

0.0 

-0.5 

4. SPAIN 

36 

8.5 

0.6 

0.0 

5. TOKYO 

35 

7.0 

-0.5 

0.4 

6. NATAL 

29 

7.4 

0.7 

0.3 

7. QUIP A 

-16 

9.5 

0.0 

-0.9 

8. SHRAZ 

29 

8.1 

-0.4 

0.7 

9. CURAC 

12 

8.1 

0.2 

-0.3 

10. JUPTR 

27 

7.4 

-0.2 

0.6 

11. VILDO 

-32 

7.6 

0.4 

0.5 

12. MAUIO 

21 

8.5 

-0.4 

-0.2 

13. AUSBK 

-31 

7.5 

0.0 

0.1 

Mean shift 


8.1 

0.0 

0.1 

rms (about mean) 

0.7 

— 

- 


Method B (Idealized Station) 


Zenith Limit 

Ah(m) 

AE(m) 

A N(m) 


7.2 

0 

0 


8.1 

0 

0 


9.3 

0 

0 
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V. SUMMARY 


The effect of omission of parallactic refraction correction for optical ob- 
servations in the determination of station coordinates was shown to have a sig- 
nificant effect on station height. Error analysis was made on a large geodetic 
satellite data distribution on 23 satellites for each of 13 Baker-Nunn camera 
sites. An average height adjustment for the 13 stations amounted to +8.1 
meters with an rms dispersion of ±0.7 meters for individual sites. The effect 
on horizontal station components for each site was less than a meter. The 
latter effect results from a good geometrical distribution of data at a station. 
Similar results for station adjustments were obtained from a satellite probability 
distribution which was assumed to represent the theoretical population for the 
above discrete satellite sample. 
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APPENDIX 

A SPECIAL CASE OF PARALLACTIC REFRACTION ERROR 


1. Modified Method A — Combined Parallactic Refraction Error 


A special case of interest is analyzed because of a program modeling error 
of parallactic refraction correction that was made in processing the Baker-Nunn 
observations for the preliminary GEM 1 solution. This case is analyzed with 
the error analysis model for Method A given in the report. The error is associ- 
ated with the projections for the parallactic angle q, given in equation (13) 
through (16) and shown in Figure 4 of the report. This material of Method A 
should be referred to for the following discussion. 

A modified parallactic angle q', defined below, was used instead of q in 
equations (15) and (16) and was applied in the opposite sense of a correction. 

This resulted in a combined error on the observation components, consisting of 
the normal refraction error as given by equations (13) and (14) plus the effect 
of the modified correction, namely 


cos 8 (c^ - a) - AR sin q + AR sin q' 


(A-l) 


§ 0 ' - 8 = - A R cos q - AR cos q' (A-2) 


The modified parallactic angle q' is similar to q which is shown in Figure 4 
of the report. The modified angle resulted by using the geocentric direction of 
the satellite in place of its topocentric direction (a, 8 ). The geometry for this 
case may be represented in Figure 4 by making the following replacements: 


replace 8 by cp s 
replace a by k s 
replace t by t 1 
replace z by 6 


geocentric latitude of satellite 

longitude of satellite 

t* = k s - k , \ station longitude 

angle between the local vertical direction of the sta- 
tion and the geocentric direction of the satellite. 


From similiarity to equations (15) and (16) in the report, 


Preceding page blank 


19 



sin q' = sin t ' cos 0/s in 9 


(A- 3) 


cos q' = (sin0 - cos 9 sin 0 s)/sin 9 sin0 s. (A-4) 


The modified observation error components, given by equations (A-l) and 
(A-2), are substituted in the matrix relation (22) of the report. A least squares 
adjustment is made for station coordinates as in the case of the normal paral- 
lactic refraction error of Method A. Results for the modified case are presented 
in the next section. 

2. Results of Modified Method A 


Table A1 presents the station adjustments due to the combined parallactic 
refraction error. The average station height adjustment is shown to be +13.5 
meters with little adjustment for the horizontal station components. A system- 
atic effect of the height adjustment A h for individual stations may be seen in 
Figure A-l where A h is plotted as a function of station latitude. This effect on 
station height corresponds to (a) nearly twice the error of that obtained for the 
normal parallactic refraction error (Table 4 of the report) for stations near the 
equator and (b) diminishes to approximately 1-1/2 times the normal case for 
stations remote from the equator. The maximum deviation in station latitude 
as shown in Table A1 occurs for SPAIN at 36° North Latitude. 

The reason for the latitude variation in A h may be seen from the difference 
between q and q' as stations deviate from the equator. The effect of the obser- 
vation errors in equations (A-3) and (A-4) on the zenith component is 


AZ = AR [1 + cos (q-q')] 


(A- 5) 


For the normal case of parallactic refraction error, AZ = AR. In the present 
case for the combined error of equation (A-5), if q is near q’ then AZ - 2 AR 
and the height adjustment is nearly twice the effect for the normal case. As q 
separates from q’ then A Z diminishes and in turn A h reduces. A qualitative 
discussion for the behavior between q and q* is given. 

For a station on the equator it may be shown, with the aid of Figure 4 and 
the associated equations for q and q», that q = q» for satellite azimuths in the 
North/South and East/West directions and that cos (q - q’) ^ 1 for values of 
azimuth in between. However, for a station at +30° latitude and for large satellite 
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Table A-l. Station Adjustments for the Modified Case of Method A 


Station 

Name 

o 

0 

A h(m) 

A E(m) 

A N(m) 

A h*(m) 

(Latitude) 

(Height) 

(Longitude) 

(Latitude) 

(Height) 

1. ORGAN 

32 

12.3 

0.4 

-1.4 

14 

2. OLFAN 

-26 

14.6 

1.0 

0.9 

7 

3. WOOMR 

-31 

16.9 

0.6 

1.6 

12 

4. SPAIN 

36 

11.5 

0.9 

-1.2 

15 

5. TOKYO 

35 

10.1 

0.0 

-0.8 

11 

6. NATAL 

29 

13.0 

-0.4 

-0.6 

15 

7. QUIPA 

-16 

16.3 

-0.4 

1.3 

7 

8. SHRAZ 

29 

12.1 

0.7 

-1.2 

15 

9. CURAC 

12 

15.3 

-1.3 

-0.4 

11 

10. JUPTR 

27 

12.5 

-1.3 

-0.3 

14 

11. VILDO 

-32 

12.7 

-0.1 

1.0 

7 

12. MAUIO 

21 

15.5 

-0.3 

-0.8 

14 

13. AUSBK 

-31 

13.9 

0.2 

0.9 

7 

Mean shift 


13.5 

0.4 

-0.8 

11.5 


* Difference between the GEM1 solution and a later solution which employed the proper 
parallatic refraction correction. 

zenith and declination angles, q and q’ separate significantly. For instance if 
azimuth is true North, then q' = 0 and q = 180° for zenith angles greater than 
60° (60° < S < 90). A similar case exists for a station at -30° latitude. The 
reason for the separation of q and q' is due to the convergent meridians toward 
the poles at S = ±90°, which affects the orientation on q as seen by Figure 4. In 
a limited region about the poles, corresponding to satellite declination near 
±90°, values of q and q' separate widely and hence AZ reduces by equation A-5. 
Since the scheduling opportunity is greater in this region for the stations that 
are closer to the poles, the reduction in A Z and thus A h is expected to be greater 
for stations that depart more from the equator as shown in Figure A-l . 

A final result is given in the last column of Table A-l. It represents the 
height differences between the preliminary GEM 1 solution and a new solution 
in which the proper parallactic refraction correction was applied. The average 
height difference for the 13 Baker-Nunn stations is 11.5 meters. This result 
may be compared with the error analysis result of 13.5 meters for combined 
parallactic refraction correction. The two mean height errors are not com- 
pletely comparable since much additional data was used in the new solution. 
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LATITUDE 


Figure A-l. Variation of A h with Station Latitude for Modified Method A 


In the case of the preliminary GEM 1 solution a nominal correction of -15 
meters in height was applied to each of the Baker-Nunn stations at the time of 
the solution. The nominal value at that time was based upon the technique of 
Method B for normal parallactic refraction error. A factor of slightly less 
than two was applied to the derived mean height adjustment from Method B in 
order to account for the effect of the combined parallactic refraction error given 
by equation A-5. 
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